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The periodically kicked quantum rotopi is 
known for non-classical effects such as quantum 
localisation in angular momentum spaced— or 
quantum resonances in rotational excitation^i^. 
These phenomena have been stndied in diverse 
systems mimicking the kicked rotor, such as cold 
atoms in optical lattices^, or conpled photonic 
strnctnres2ri£. Recently, it was predictedii — 
that several solid state quantnm localisation phe¬ 
nomena Anderson localisation^^, Bloch oscilla- 
tion a^^d^ , and Tamm-Shockley surface statesiid^ 
— may manifest themselves in the rotational dy¬ 
namics of laser-kicked molecules. Here, we report 
the first observation of rotational Bloch oscilla¬ 
tions in a gas of nitrogen molecnles kicked by a 
periodic train of femtosecond laser pulses. A con¬ 
trollable detuning from the quantum resonance 
creates an effective accelerating potential in an¬ 
gular momentum space, inducing Bloch-like oscil¬ 
lations of the rotational excitation. These oscilla¬ 
tions are measured via the temporal modulation 
of the refractive index of the gas. Onr resnlts in¬ 
troduce room-temperature laser-kicked molecules 
as a new laboratory for studies of localisation phe¬ 
nomena in quantum transport. 

Bloch oscillations present one of the most famous 
and intriguing quantum localisation effects in solid state 
physics: Electrons in crystalline solids subject to an ex¬ 
ternal dc electric field exhibit oscillatory motion instead 
of a mere uniform acceleration as in empty space. This 
effect was predicted at the inception of quantum mechan¬ 
ics in 192 9^^d^ , but it took more than 60 years to first 
observe Bloch oscillations in semiconductor superlattice 
structures'^. Bloch oscillations were later observed in 
the momentum distribution of ultracold atoms driven 
by a constant force^S. Recently, it was proposed to in¬ 
duce Bloch oscillations in the rotation of linear molecules 
excited by periodic trains of short nonresonant laser 
pulsesi^. This proposal is based on an analogy between 
coherent rotational excitation caused by the laser pulses 
and propagation of an electron in a one-dimensional pe¬ 
riodic lattice. 

A non-resonant, linearly polarised short laser pulse 
acts as a kick to a molecule. It excites a rotational 
wave packet via multiple coherent Raman-type inter- 
actions2ii2^: |d>)(t) = Cj{t)\J), where |J) are the 

angular momentum states, and the projection quantum 
number M is dropped since it is not changing. Under 
field-free conditions, the coefficients Cj{t) oscillate as 
exp[—iBJ{J + where B is the molecular rota¬ 


tional constant. The dynamics of the wave packet are 
determined by a single parameter, the rotational revival 
time Aev = Trh/B. The wave packet revives exactly at 
integer multiples of the revival time^^ when all the time- 
dependent phase factors become equal to unity. 

Consider a train of laser kicks with a constant time 
delay r between the pulses. The time delay is chosen 
to be slightly detuned from the rotational revival time: 
T = (1 -I- J)trev Due to the detuning, each compo¬ 
nent of the rotational wave packet acquires a small phase 
4>j = TrSJ{J + 1) from pulse to pulse (integer multiples 
of 27r are dropped). From now on, we follow the dy¬ 
namics stroboscopically, by considering only Cjin), the 
wave packet coefficients just after the n**' kick. For weak 
pulses and small detuning, the change of Cj{n) from the 
to the (n -I- 1)*^ pulse is given as (see supplementary 
material) 

Cj{n+l)-Cj{n) ^ 

ij [Cj+ 2 in) + Cj. 2 {n)] - icj>{J) Cj{n ). (1) 

The first term on the right hand side describes the 
laser coupling of the rotational levels, where P = 
(Aa/Ah) f E'^{t)dt is an effective strength of the laser 
pulse^i (Aa is the molecular polarisability anisotropy, 
and E{t) is the envelope of the laser electric field). Note 
that the laser field only couples states of the same parity, 
AJ = 0, ±2. The second term is caused by the detuning 
of the pulse train period from the rotational revival time 
Aev As the change of Cj{n) is small, the difference equa¬ 
tion © can be recast as a differential equation, where n 
becomes a continuous variable: 

= -^[Cj+ 2 {n) + Cj. 2 {n)] +</.(J)Cj(n). (2) 

Equation ([2]) looks like the Schrodinger equation for a 
particle moving in a periodic ID lattice, where the num¬ 
ber of pulses n plays the role of dimensionless time. The 
sites of the lattice are the angular momentum states | J). 
(Note that even and odd J form two independent lat¬ 
tices.) The first term on the right hand side of equa¬ 
tion @ describes the coherent hopping between lattice 
sites. The second term can be interpreted as an on-site 
potential V{J) = 4>{J) = tt5J{J + 1). 

At the quantum resonance, t = Uev (i-e. 5 = 0), the 
potential is V (J) = 0. In this case, equation ([2]) describes 
the standard tight-binding model in solid state physics^^, 
depicted in figure [1] (a). The eigenstates of this model 
are Bloch waves, which are characterised by the continu¬ 
ous quasimomentum k and the energy dispersion relation 
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a Quantum Resonance: tight-binding modei 
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FIG. 1. The periodically kicked rotor as a particle 
in a ID lattice, a, For a rotor kicked at the exact quan¬ 
tum resonance, the dynamics can be described by a tight- 
binding model. The angular momentum levels J form a dis¬ 
crete ID grid, and the laser pulses couple sites with A J = ±2, 
where the coupling strength is proportional to the effective 
kick strength P of the pulses (see text), b, A detuning from 
the quantum resonance introduces an effective potential V{J) 
to the model. The dynamics are then similar to Bloch oscilla¬ 
tions: A wave packet (particle) is accelerated by the effective 
potential, and its wave vector k grows. When it reaches the 
edge of the Brillouin zone (here at k = 7r/2), the wave is 
reflected by Bragg reflection. 


e(fc). The latter is given as £{k) = —(P/2) cos(2/c). Wave 
packets formed by these states propagate without limita¬ 
tion, reaching very large J values. This is the signature 
of the rotational quantum resonance effect^i^. 

To treat the case of non-zero detuning, one can derive 
semi-classical equations of motion for the wave vector k 
and the lattice coordinate J— : 


dk dy(J) dJ de(fc) 


dn 


dJ dn 


dk 


= Psin(2fc). (3) 


The first equation is Newton’s second law, and the sec¬ 
ond one defines the group velocity of the Bloch waves in 
the lattice. Equations (| 31 ) can also be derived in a non- 
perturbative way for strong laser pulses using the formal¬ 
ism of e-classics2^. For negative detuning 5 < 0, V{J) is 
an accelerating potential [as depicted in figured] (b)], and 
equations (I3|) are similar to the ones describing electrons 
in crystalline solids subject to a constant electric field, the 
famous problem treated by Bloch and Zener— For low 
J, the force —dV/dJ is weak, and the dynamics resembles 
the case of the quantum resonance. As the Bloch wave 
packet moves to larger J, the quasimomentum k grows 
and approaches the end of the Brillouin zone at fc = 7r/2, 
where the length of the Bloch waves is comparable to 
the “lattice constant” AJ = 2. As a result, the wave is 
reflected due to Bragg scattering. The now backwards 
moving wave packet is decelerated by the potential, until 
it stops and the cycle starts again. 

To illustrate the above predictions, we simulated the 
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FIG. 2. Bloch oscillations in the angular momentum 
of laser kicked ^‘*N 2 molecules. Simulated population of 
the angular momentum levels J, as a function of the number 
of laser pulses. The initial rotational temperature is 298 K. 
The considered pulse train parameters are r = 8.36 ps (0.2% 
less than the revival time of 8.38 ps) and an effective interac¬ 
tion strength of P = 5 (see text). 


population distribution of the rotational levels J for 
molecules at room temperature, interacting with 
a train of pulses with P = 5. The detuning is equal to 
S = —0.2%. The result of the simulation is shown in fig¬ 
ure [2| One can clearly see the predicted oscillations of 
the rotational excitation. For the first four pulses, the 
angular momentum population shifts to higher J; the 
shift per pulse is constant. From the the fifth pulse on, 
the movement is reversed and directed towards lower J. 
After eight pulses, the system returns approximately to 
the initial state, and the cycle starts again. 

To observe the rotational Bloch oscillations experimen¬ 
tally, we employ a scheme similar to the one proposed 
earlier—. A triple nested interferometer set-up is used to 
create a periodic train of eight femtosecond laser pulses^i. 
These pulses interact with nitrogen molecules at room 
temperature and standard pressure conditions. The ro¬ 
tational excitation of the molecules leads to molecular 
alignment, which causes optical birefringence. We mea¬ 
sure the time-dependent alignment signal and take its 
average over one rotational revival time as a measure of 
the rotational excitation (see Methods). 

Figure [3| shows the results for several pulse train pa¬ 
rameters P and (5, and demonstrates the predicted oscil¬ 
latory behaviour. The following trends can be observed: 
First, the amplitude of the oscillations increases with the 
effective kick strength P, but decreases when the detun¬ 
ing |i5| is increased. In addition, the oscillation period de¬ 
creases both with increasing P and |<5|. The same trends 
are found for the solutions of equations ([3]). The statis¬ 
tical uncertainties are smaller than the plot symbols. 

To further compare the predictions with the experi¬ 
ments, we solved equations choosing J(0) =0 and 
k{0) = 7r/4 as initial conditions. This value of the 
quasimomentum corresponds to an initial growth rate 
of dJ/dn = P of the angular momentum (the typical 
change of J by a single pulse). Figure Sj shows the oscil- 
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FIG. 3. Measured and simulated time-averaged align¬ 
ment signal. A comparison of the measured and simulated 
time-averaged alignment signals for periodically kicked nitro¬ 
gen molecules, for different values of the kick strength P and 
the detuning 5. The statistical error of the experimental data 
is much less than the plot symbols and therefore not plotted. 
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FIG. 4. Period of the Bloch oscillation. The oscillation 
period as a function of the detuning, for different effective 
interaction strengths P. The lines are the solution of the semi- 
classical model, equations ®, the markers are experimental 
values. 

lation period as a function of the detuning 5 for different 
kick strengths. The curve predicted by the semi-classical 
model (solid lines) and the values extracted from the ex¬ 
periment (markers) show general agreement. The visibil¬ 
ity of the measured oscillations is reduced compared to 
the simulations presented in figure |31 We attribute this 


to collisional decoherence as well as to imperfections of 
the pulse train. 

Summarising, we described the first observation of 
Bloch oscillations in a quantum rotational system. We 
used nitrogen molecules at room temperature and stan¬ 
dard pressure conditions, kicked by a periodic train of 
femtosecond laser pulses, close to the quantum resonance 
condition. Our results demonstrate that laser kicked 
molecules are subject to the same localisation phenomena 
that are seen in quantum transport. 


METHODS 

Experiments. The data are measured using a stan¬ 
dard ultrafast pump probe set-up, depicted in figure [S] 
A triple-nested interferometer set-up is utilised to split 
a single 3 mJ 70 fs 800 nm pulse into an periodic eight 
pulse pump train, as well as a weak probe pulse as de¬ 
scribed previously^!. This train is focused by an off-axis 
parabola into a fused silica cell containing a constant flow 
of dry gas. The induced rotational dynamics are en¬ 
coded in the time-dependence of the angular alignment 
of the molecules (for recent reviews on laser molecular 
alignment see!^“— ), and gives rise to a weak optical bire¬ 
fringence. This time-dependent birefringence is probed 
using a circularly polarized 400 nm probe pulse with a 
variable optical path length. The birefringence within 
the sample causes the probe’s polarisation to become 
weakly elliptical after propagation. The probe is split 
with a Glan-type polariser, and the polarisation compo¬ 
nent amplitudes are recorded as a function of probe de¬ 
lay time. This measurement is performed with the pump 
train both present and blocked in rapid succession, and 
the normalised difference is proportional to the molecu¬ 
lar alignment signal at the selected time delay between 
the pump and probe arms. A detailed description of the 
set-up is presented in precedent work!^. 

We use the fact that the time-averaged alignment is 
effectively a monotonic function of the angular momen¬ 
tum J, especially at high temperatures!^. An increase 
(decrease) in the expectation value (J) of the angular 
momentum translates into an increase (decrease) of the 
time-averaged alignment. Therefore, one can observe the 
Bloch oscillations indirectly via the time-averaged align¬ 
ment signal (population alignment). Each value of the 
population alignment is found by averaging the molecu¬ 
lar alignment signal over one rotational revival time; this 
signal is sampled at about 100 uniformly spaced probe 
delay times. The probe times are arranged to accurately 
show the population alignment and avoid contamination 
by the coherent alignment, which has different decoher¬ 
ence mechanisms. Each population alignment value is 
measured over 150, 000 laser shots and has negligible sta¬ 
tistical error. The oscillation periods presented in figure^] 
are obtained from a polynomial least square fit of the 
time-averaged alignment data. 

Simulations. We model the non-resonant laser- 
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FIG. 5. Simplified sketch of the experimental set-up. A Ti:Sapph source generates an 800 nm pulse that enters a 
pump-probe set-up. The pump is split into an eight pulse train using three nested interferometers. The probe pulse is converted 
to circularly polarised 400 nm with variable path length. The beams are focused using an off-axis parabola (depicted as a lens). 
Molecular alignment triggered by the strong pump train causes birefringence that alter the circular polarisation of the weak 
probe pulse. The probe’s polarisation is split and measured, yielding a time-dependent molecular alignment signal that can be 
time-averaged to yield population alignment. The abbreviations used are: BS - beam Splitter, IF - interferometer, HWP - half 
wave plate, L - lens, DL - automated delay line, BBO - /1-BBO crystal, QWP - quarter wave plate, BPF - blue pass filter, GC 
- gas cell, GTP - Gian-Taylor polariser, PD - photodiode. 


molecule interaction by the effective potential U = 
— (Aa/d) cos^ where Aa is the molecular polar- 

isability anisotropy, and E[t) is the envelope of the laser 
electric held, and 9 is the angle between the molecular 
axis and the laser polarisation direction. The wave func¬ 
tion is expanded in the spherical harmonics, and we solve 
numerically the time-dependent Schrodinger equation to 
obtain the expansion coefficients. To take into account 
thermal effects, we average over the initial states, where 
each result is weighted by the Boltzmann factor (includ¬ 
ing nuclear spin statistics) of the initial state. A de¬ 
tailed description of the numerical method can be found 
in Ref.-^. 
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